Simulations of atomic trajectories near a dielectric 
surface 



N. P. Stern^ D. J. Alton^ and H. J. Kimble^ 

^ Norman Bridge Laboratory of Physics 12-33, California Institute of Technology, 
Pasadena, California 91125, USA 

E-mail: hjkinible@caltech.edu 

Abstract. We present a semiclassical model of an atom moving in the evanescent 
field of a microtoroidal resonator. Atoms falling through whispering-gallery modes can 
achieve strong, coherent couphng with the cavity at distances of approximately 100 
nanometers from the surface; in this regime, surface-induced Casmir-Polder level shifts 
become significant for atomic motion and detection. Atomic transit events detected 
in recent experiments are analyzed with our simulation, which is extended to consider 
atom trapping in the evanescent field of a microtoroid. 



PACS numbers: 37.30.-|-i , 34.35.-ha, 42.50.Ct, 37.10.Vz 



Simulations of atomic trajectories near a dielectric surface 



2 



1. Introduction 

Strong, coherent interactions between atoms and light are an attractive resource for 
storing, manipulating, and retrieving quantum information in a quantum network with 
atoms serving as nodes for quantum processing and storage and with photons acting as 
a long-distance carrier for communication of quantum information [1]. One realization 
of a quantum node is an optical cavity, where light-matter interactions are enhanced 
by confining optical fields to small mode volumes. In the canonical implementation, 
a Fabry-Perot resonator with intracavity trapped atoms enables a panoply of cavity 
quantum electrodynamics (cQED) phenomena using single photons and single atoms, 
and thereby, validates many aspects of a cQED quantum node [H |3] . 

Despite these achievements, high-quality Fabry-Perot mirror cavities typically 
require significant care to construct and complex experimental instrumentation to 
stabihze. These practical issues have begun to be addressed by atom chips [11[5], in which 
atoms are manipulated in integrated on-chip microcavity structures offering a scalable 
interface between light and matter [6l[7l|8]. Owing to their high quality factors, low mode 
volumes, and efficient coupling to tapered optical fibers |D] , microtoroidal resonators are 
a promising example of microcavities well-suited for on-chip cQED with single atoms 
and single photons [10] • Strong coupling [HI [12] and non-classical regulation of optical 
fields [131 E] have been demonstrated with atoms and the whispering-gallery modes of 
a silica microtoroidal resonator. 

In our experiments with microtoroids, Cs atoms are released from an optical trap 
and fall near a silica toroid, undergoing coherent interactions with cavity modes as each 
atom individually transits through the evanescent field of the resonator. In the most 
recent work of [12j, atom transits are triggered in real-time to enable measurement 
of the Rabi-split spectrum of a strongly-coupled cQED system. Whereas a single 
atom is sufficient to modify the cavity dynamics, falling atoms are coupled to the 
cavity for only a few microseconds. Atom dropping experiments necessarily involve 
a large ensemble of individual atomic trajectories and represent, consequently, a far 
more complex measurement result. 

Interactions between a neutral atom and a dielectric surface modify the radiative 
environment of the atom resulting in an enhanced decay rate [I5] and Casimir-Polder 
(CP) forces [161 EZ]- These perturbative radiative surface interactions are usually 
insignificant in cQED experiments with Fabry-Perot resonators where atoms are far from 
mirror surfaces, but in microcavity cQED, atoms are localized in evanescent fields with 
scale lengths A/27r ~ 150 nm near a dielectric surface. The experimental conditions for 
microtoroidal cQED with falling atoms in [T2| necessarily involve significant CP forces 
and level shifts while simultaneously addressing strong coupling to optical cavity modes. 
Theoretical analysis of this experiment requires addressing both the strong atom-cavity 
interactions and atom interactions with the dielectric surface of the microtoroid. As 
reported in [12], spectral and temporal measurements offer signatures of both strong 
coupling to the cavity mode and the significant infiuence of surface interactions on 
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atomic motion. The role of these effects is quantified with detailed simulation of the 
trajectories of falling atoms detected in the real-time at low photon numbers. 

In this article, we discuss in detail the approach used to simulate atomic motion 
near the surface of an axisymmetric dielectric resonator under the influence of strong 
coherent interactions with cavity modes. The experimental detection method of |T2] 
is implemented stochastically in a semiclassical simulation of atom trajectories. These 
simulations provide a perspective on the atomic motion of atom transits recorded in our 
microtoroid experiments, while offering additional insights into the loading of optical 
evanescent field traps. In section [2| we outline the semiclassical model of a two-level 
atom coupled to the whispering gallery modes of a microtoroidal resonator. In section [3} 
we review the optical dipole forces which are a critical factor influencing atomic motion 
in an optical cavity. Our calculations of modified emission rates and Casimir-Polder 
surface interactions are detailed in section |4| Section [5] describes the implementation 
of our model for simulating recent atom-toroid experiments. Finally, section [6] extends 
our simulation to evanescent field traps around a microtoroid. 

2. Atoms in a microtoroidal cavity 

We approach the motion of atoms moving under the influence of surface interactions 
and coherent cavity dynamics with a semiclassical method to efficiently simulate a large 
number of atom trajectories. For surface interactions, dispersion forces are calculated 
perturbatively using the linear response functions of Si02 and a multi-level atom. For 
nearly-resonant non-perturbative coherent interactions between atom and cavity, the 
atomic internal state and the cavity field are treated quantum mechanically within the 
two-level and rotating-wave approximations. 

Simulations of atomic motion follow the semiclassical method detailed in [18]. 
Mechanical effects of light are incorporated classically as a force F{r) on a point particle 
atom at location r. Trajectories f{t) are calculated with a Langevin equation approach 
to incorporate momentum diffusion from fluctuations. At each simulation time step t*, 
the atomic velocity is calculated as: 



where is the velocity at the i time step, mcs is the atomic mass, and At is the 



standard deviation of 1. Given the force F and diffusion tensor Dij as discussed in 
section |3| the atom trajectory r(t) and cavity transmission and reflection coefficients, 
T(t) and R(t) are calculated. A single atom strongly coupled to the cavity mode has a 
large effect on cavity fields and optical forces, requiring simultaneous solutions of atomic 
motion and cQED dynamics. 

Full quantization of atomic motion leads to an unwieldy Hilbert space not conducive 
to efficient simulation. In contrast, semiclassical methods are well-suited for simulating 
atomic motion in experiments with falling atoms near resonators. The ratio of the recoil 
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simulation time step f 



— f. The VT' are normally distributed with zero mean and 
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energy to the linewidth of the cesium 65*1/2 QP3/2 transition is less than 10~^. Further, 
the recoil velocity of ~ 3.5 mm/ s is much less than the typical velocity of falling atoms of 
order 200 mm/s so that each spontaneous emission event represents a small momentum 
kick. Cavity fields and internal atomic states respond quickly to environment changes, 
allowing calculation of optical forces and momentum diffusion in a constant-velocity 
limit at time t and energy shifts from surface interactions as if atom the atom were 
stationary. The remainder of this section discusses the quantum mechanical equations 
of motion for the atom and cavity fields in the low-probe intensity limit, to be followed 
later by contributions to the force F used in ([l). 

2.1. Modes of a microtoroidal resonator 

An idealized microtoroid has axial symmetry, so we work in a standard cylindrical 
coordinate system r — )■ (p, </), z). The toroid is modeled as a circle of diameter Dm 
with dielectric constant e revolved around the 2;-axis to make a torus of major diameter 
Dm (Fig- The toroid is therefore defined by its minor diameter and its 

principal diameter Dp = Du + D^. The fabrication and characterization of high-quality 
microtoroids are described in detail elsewhere [9]. 

The axisymmetric cavity modes of interest are whispering-gallery modes (WGM), 
which lie near the edge of the resonator surface and circulate in either a clockwise 
or counter-clockwise direction. These modes are characterized by an azimuthal mode 
number m, whose magnitude gives the periodicity around the toroid and whose 
sign indicates the direction of propagation. The WGMs for ±m are degenerate in 
frequency but travel around the toroid in opposite directions. The mode electric 
fields for the WGM traveling waves are written as E{f) = E^s_xf{p,z)e^"^'^, where 
f{p,z) = E{p, z)/ E-^aax is the mode function in the p — z cross-section normalized by 
the maximum electric field -Emax- In general, backscattering couples these two modes 
so that a more useful eigenbasis for the system consists of the normal, standing wave 
modes characterized by a phase and the periodicity \m\. This backscattering coupling 
h is assumed to be real, with the phase absorbed into the definition of the origin of 
the coordinate 0. In addition, the mode's field decays at a rate through radiation, 
scattering, and absorption. In our simulations, a cavity mode is fully described by its 
spatial mode function /(r), its azimuthal mode number m, its loss rate n^, and the 
coupling h to the counter-propagating mode with mode number — m. 

We model the microtoroid modes using a commercial finite-element software 
package (COMSOL) to solve numerically for the vector mode functions /(p, z) for modes 
of a given m [19]. Mode volumes are calculated from. 



JdVe{f)\E{f^^^ 

max 

In this notation [10], the coupling of a circularly polarized optical field to an atomic 



Vm = ^ ^ ^' =2n ldAe{p,z)pf{p,zy (2) 
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dipole located in the evanescent field of the cavity is calculated as: 

gir-) = {d- E) = fip,z)e-\[^^ (3) 

where cTis the dipole operator and a;i°^ = 2ttc/Xo is the vacuum transition frequency of 
the two- level atom with free-space wavelength Aq. WGMs are predominantly linearly 
polarized, and so we average over the dipole matrix elements to obtain an effective 
traveling wave coupling gtw which is approximately ~ 0.6 of the value for circularly 
polarized light (see supplementary information of [11]. Travelling wave modes of an 
axisymmetric resonator are not strictly transverse. For the toroid geometries considered 
here, with Dp, ^ A, the azimuthal component is small and we assume that the 
optical field is linear outside of the toroid. Since the cavity losses are dominated by 
absorption and defect scattering rather than the radiative lifetime set by the toroid 
geometry pTO] , we let Ki and h be experimental parameters. Fig. [l] shows the lowest- 
order mode with m = 118 for a toroid with {Dp,Dyyi} = {24,3} fim. The index m is 
chosen so that the cavity frequency oUc is near the 65*1/2 QP3/2 transition of Cs with 
J^^/2n = 351.7 THz. 

The local polarization of modes varies throughout the interior and exterior of the 
toroid. Approximate solutions for constant polarization suggest classifications as quasi- 
transverse modes, labeled transverse electric (TE) and transverse magnetic (TM) modes, 
although actual solutions are not transverse. A reasonable analytic approximation for 
the lowest-order mode function with mode number m outside of the toroid is that of a 
Gaussian wrapped around the toroid's surface that decays exponentially with distance 
scale set by the free space wavevector I/Aq = 27r/Ao, 

/(p,^) ~ e-'^/^°e-('^/'^o)' (4) 

where d{p, z) = \/{p — -Dm/2)^ + z"^ — -Dm/2 is the distance to the toroid surface, 
ip{p,z) = arctan p_^^/2 angle around the toroid cross-section {ip = at z = 0), 

and ipQ is a characteristic mode width (see Fig. la). Higher order angular modes are 
characterized by additional nodes along the coordinate ip- 



2.2. Cavity QED in an axisymmetric resonator 

We consider a quantum model of a two-level atom at position r{t) coupled to an 
axisymmetric resonator shown schematically in Fig. |2} The terminology used here 
follows the supplemental material of [H], [13], and [12], but the general formalism 



can be found in additional sources (see [20], for example). As described in section 2.1 
an axisymmetric resonator supports two degenerate counter-propagating whispering- 
gallery modes at resonance frequency Uc to which we associate the annihilation (creation) 
operators a and h {^a) and Ij^). Each traveling-wave mode has an intrinsic loss rate, 
the modes are coupled via scattering at rate h. External optical access to the cavity is 
provided by a tapered fiber carrying input fields {om, fem} at probe frequency Up. Fiber 
fields couple to the cavity modes with an external coupling rate Kex- The output fields of 




m 



Figure 1. (a) A scanning electron microscope image of a microtoroid with definitions 
of the relevant parameters discussed in the text, (b) The lowest order mode function 
z) of a toroid mode with {Up, i^m} = {24, 3} ^vci and m = 118 and A = 852 mn. 



the fiber taper in each direction are the coherent sum of the input field and the leaking 
cavity field, {oout, &out} = -{oin, hm] + y2fi~{a, 6} [m Il3]. 

We specialize to the situation of single-sided excitation, where (6in) = 0. The 
input field Oin drives the a mode with strength = iy^2K^{ain) so that the incident 
photon flux is Pin = (ajj^ajn) = |epp/2/€ex- Experimentally accessible quantities are 
the transmitted and reflected photon fluxes, Pt = ('^Lt'^out) and Pr = (^lut^out), 
respectively. In experiments, data is typically presented as normalized transmission 
and reflection coefficients defined as T = Px/Pin and R = P-q^/Pi^. In the absence of an 
atom, the functions T and R for the bare cavity depend on the detuning Acp = — Wp 
and the cavity rates /i, Ki, and k^^. At critical coupling, /tex = a/k^ + h?., the bare cavity 
T ^ when Aep = [21]. 

The cavity modes {a, 6} both couple to a single two-level atom with transition 
frequency u^, at location f. In the context of cQED, the atomic system is described by 
a single transition with frequency cua with the associated raising and lowering operators 
cr^ and a~ and an excited state field decay rate 7. The atomic frequency ujs_{r) may 
be shifted from the free-space value wi^^ by frequency 5a (^) due to interactions with the 
dielectric surface. The coupling of the traveling-wave modes {a, 6} to the atomic dipole 
is described by the single-photon coupling rate gt^nif) = Qt^"^ fi.Pi ^)^^^^ 1 where f{p,z) 
is the cavity mode function and 9 = mcj). A discussion of /(p, z) for the modes of a 
microtoroid appears in section 2.1, For an atom in motion, uji^{r), 7(r), and gty^{r) are 
spatially-dependent quantities that depend on the atomic position r{t). 

To study the atom-cavity dynamics, we write the standard Jaynes-Cummings-style 
cQED Hamiltonian for coupled field modes [221 [TT] : 

H/h = u^{f)(x+a- + {a^a + b^b) + h {a^b + b^a) + {ele'^^^'a + Spe-'^'^'a^) 
+ {g;^ir)a^a- + (7tw(r)a+a) + ((7tw(r)&V- + gUf^cr+b) . (5) 
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Figure 2. (a) Schematic of the atom-toroid system. Coherent optical fields in the 
tapered fiber couple into whispering-gallery cavity modes of an axisymmetric resonator. 
These fields can couple to an atomic transition with rate g, scatter to the counter- 
propagating mode (h), escape to the environment (^i), or couple in/out of the fiber 
(kox)- The atom is described as a two-level system with transition frequency and 
spontaneous emission rate 7. (b) Imaginary part of the eigenvalues of the linearized 



systems as a function of detuning A = 



,(0) 



for a Cs atom at 



g = 60 MHz critically coupled to a cavity with parameters {Ki,/i}/27r 



= •7r/4 and 
{8,0} MHz 



(Eqs. (9a)) 



Following the rotating-wave approximation, we write the Hamiltonian in a frame 
rotating at [m UHl ED] : 

H/h = Aap(r)o-+o-" + Acp {a^a + b^b) + h (a^b + b^a) + e*a + SpU^ 

+ ((7L(f)aV- + (7tw(r)(T+a) + ((7tw(r)6V- + g^f^a+b) , (6) 

where Aap(r) = 0Ji,{r) — Up and Acp = — Wp. Dissipation from coupling to external 
modes is treated using the master equation for the density operator of the system p: 

p = — p] + ^ (2apa^ — a'^ap — pa'^a^ + k (26p6^ — 6^6p — p6^6) 

+ 'J (2a^ pa^ — p — pa^ a^) (7) 

Here, k = Ki + Kex is the total field decay rate of each cavity mode, and 2'y{f) is the 
atomic dipole spontaneous emission rate, which is orientation dependent near a dielectric 
surface (section 4.1). 

The Hamiltonian ^ can be rewritten in a standing wave basis using normal modes 
A= (a + b) /V2 and B = {a-b) /v^, 

H/h = Aap(f)a+a- + (Aep + h)A^ A + (A^p - h)B^B 
+ [elA + EpA^) /V2 + (45 + EpB^) /V2 

+ gp,{r) [X^a- + a+ A) - ig^{r) [B^a' - a+B) , (8) 

where gAir) = g^^^f (p, z) cos 9, g^ir) = g^s^^f (p, z) sin 9, and ^f^ax = V^gZ^""- In 
the absence of atomic couphng {gtv, = 0), these normal modes are eigenstates of 
With gtvj 7^ 0, the eigenstates of the Hamiltonian are dressed states of atom-cavity 
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excitations. With h = and gtw 7^ 0, the atom defines a natural basis in which it 
couples to only a single standing wave mode. For the modes {A,B} defined above, 
coupling may occur predominantly, or even exclusively, to one of the two normal modes 
depending the azimuthal coordinate 6. For such 6, the system can be interpreted as an 
atom coupled to one normal mode in a traditional Jaynes-Cummings model with dressed- 
state splitting given by the single-photon Rabi frequency = 2g = 2g^a.xf{p, z), along 
with a second complementary cavity mode uncoupled to the atom. Approximately for 
gtw ^ h, this interpretation is consistent for any arbitrary atomic coordinate 6. For 
h ^ and comparable to Hi with a fixed phase convention (such as lm.{h) = used 
here), this decomposition is not possible for arbitrary atomic coordinate 6; the atom in 
general couples to both normal modes as a function of (p pT| [20] . 

The master equation ([T]) can be numerically solved using a truncated number 
state basis for the cavity modes. Alternatively, the system is linearized by treating 
the atom operators as approximate bosonic harmonic oscillator operators with 
[cr^,cr+] ^ 1. For a sufficiently weak probe field, the atomic excited state population 
is small enough that the oscillator has negligible population above the first excited 
level and the harmonic approximation is quite good. As part of this linearization, 
we factor expectation values of normally-ordered operator products into products 
of operator expectation values [IHl 123] • Reducing operators to complex numbers 
suppresses coherence, but numerical calculations confirm that this approximation is 
accurate when calculating cavity output fields and classical forces for the weak driving 
power levels considered here. In particular, experiments typically utilize a photon fiux 
Pt = (ajjut'^out) ~ 15 cts//is corresponding to an average cavity photon population of 
(a^a) ~ 0.1. At these photon numbers, cavity expectation values effectively factorize 
such that (a'^a) ^ (a'') (a) for the semiclassical treatment used here [21]. We use this 
approximation to write Pt = ('^Lt'^out) ~ (c^Lt) ('^out), implying that we only need 
the complex number (flout) = — («in) + A/2Kex(a) and its conjugate to calculate the 
cavity transmission at these photon numbers. This approximation is not sufficient 
for calculation of the g^'^^r) correlation function where the nonlinearities must be 
included |12J. 

The relevant equations of motion for the field amplitudes of the linearized master 
equation are, 

(a) =_(«; + iAcp) (a) - ih (b) - ie^ - ig*^ (a"), (9a) 

(6) = -{k + iAcp) (&) - th (a) - tgt^ {a'), {9b) 

{(^') = -{1 + ^Aap) {a') - igtw {a) - igl^ {b). (9c) 

Time and spectral dependence of this system of equations are governed by its 
eigenvalues Aj. The imaginary part of the eigenvalues as a function of detuning 
A = Acp — Aap = cOc — coa. are illustrated in Fig. |2|d. For large A ^ \gtw\, the three 
eigenvalues include one atom-like eigenvalue and two cavity-like eigenvalues separated by 
the mode splitting h. For intermediate A, there is an anti-crossing of two dressed-state 
eigenvalues A±, while the third (cavity-hke) Aq is uncoupled to the atom. 
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For a slowly-moving atom, the mode fields remain approximately in steady state 
as the parameters evolve with the atom trajectory r{t). Analytic steady-state solutions 



to (9a) for (a)ss and (6)ss are: 



(a)ss 

(fe)ss 



le 



(7 + iAap) [(k + iAcp) (7 + iAap) + Ifftv 



\ih(ci + iAap) + (5t*w)^] {^Kl + iAap) + gL] - [('^ + ^Acp)(7 + iAap) + |5twP]^ 
i/l (7 + iAap) + 5tw 



(k + iAcp)(7 + iAap) + |5tv 

ss 



r(a>s 



7 + iA 



(10a) 

(106) 
(10 c) 



ap 



3. Optical forces on an atom in a cavity 

Neutral atoms experience forces from the interaction of the atomic dipole moment with 
the radiation field. These optical dipole forces have a quantum mechanical interpretation 
as coherent photon scattering [221 • For a light field near resonance with the atomic 
dipole transition, these optical forces can be quite strong, even at the single photon level; 
cavity-enhanced dipole forces [IHl [22] have been exploited to trap |2H] and localize [22] 
a single atom with the force generated by a single strongly-coupled photon. In this 
section, we discuss how the optical forces, their first-order velocity dependence, and 
their fiuctuations are included in our semiclassical simulation. 



3.1. Dipole forces 

In a quantum mechanical treatment of light-matter interactions the eigenstates 
of the system are dressed states of atom and optical field. The quantum mechanical 
optical force on the atom at location r can be found from the commutator of the atom 
momentum p with the interaction Hamiltonian i^int consisting of the last two terms 
from the Hamiltonian (|6]): 

F=Jt={ [^int,Pl = -W(7L(f) [a^a- + - W(7tw(r) (a+a + h^^a') (11) 

The gradient from the position space representation of the momentum operator p only 
acts on gtw{r) and not on the field operators [30l [Slj. The steady-state expectation 



values of (11) give the dipole force on the atom in the semiclassical approximation: 

(F),, = - hVglif) ((at),,(a-),s + ((T+)ss(&)ss) 

- W^7tw(r) ((^+)ss(a)ss + {b^)ss{cr-)ss) (12) 

As described in section [2} the steady-state operator expressions are simplified by 
reducing expectation values of operator products to products of linearized steady-state 
operator expectation values. Ignoring fiber and spontaneous emission losses, an effective 



conservative dipole potential Ud can be defined by integration of (12). 
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3.2. Velocity-dependent forces on an atom 



Non-zero velocity effects on the force (12) are found by including a first-order velocity 



correction in the steady state expectation values [23l |25l [3T]. Consider a vector of 



operators O whose expectation values obey a linearized equation system such as (9a). 
Assuming a small velocity, we expand the operator expectation values (O) as: 

(d) = (d)o + (d)i + ..., (13) 

where the subscripts denote the order of the velocity v in each term. If an atom is 
moving through these fields, then the cavity parameters depend in general on atomic 
position r. As r changes in time, the fields must evolve in response. Consequently, the 
time derivative of the expectation value evolves not only from explicit time dependence, 
but from atomic motion as well. 

{d)=(^^+v-vyd) (14) 



Setting the explicit time derivatives in (14) to zero, the perturbative expansion of the 



time derivative can be equated to the original linearized equation system. Collecting 
terms of each order in velocity gives an equation for the first-order term (0)i in terms 
of the zero-velocity steady-state solution (0)o- This procedure requires the spatial 
derivative of the zero-order steady-state solutions, where spatial dependence enters 
through the atomic transition frequency Wa(r), the spontaneous emission rate 7(r), and 
the atom-cavity coupling g{f). Only terms linear in velocity are kept in the operator 



products of the force F{f) in (12) 



In practice, first-order velocity corrections are small in our simulation. For example, 
Doppler shifts arising from spatial derivatives of the cavity modes are on the order of 
k ■ V, where k is the mode wavevector. For typical azimuthal velocities of less than 
0.1 m/s, the Doppler shift is less than 1 MHz. The effect becomes more significant as 
atoms accelerate to high velocities near the surface, but atomic level shifts from surface 
interactions are more significant in this regime than the Doppler shifts. 

3.3. Momentum diffusion and the diffusion tensor in a cavity 

Quantum fluctuations of optical forces are treated by adding a stochastic momentum 
diffusion contribution to the atomic velocity in the Langevin equations of motion. We 
calculate the diffusion tensor components used in ([T]), Da, using general expressions for 
diffusion in an atom-cavity system generalized for the two-mode cavity of a toroid [32] '■ 

2Du = {hky2^\{a-)J + I (a-) J 2^ + 2k (|W. {a) J + |W. {b)J) (15) 

for i = x,y, z, where 7 is the atomic field spontaneous decay rate. The first 
term represents fluctuations from spontaneous emission, the second term describes a 
fluctuating atomic dipole coupled to a cavity field, and the third represents a fluctuating 



cavity field coupled to an atomic dipole. (15) is approximated using steady-state fields 



calculated from the linearized solutions to the master equation ( 10a). Although included 
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in the trajectory model, momentum diffusion does not significantly alter averages over 
ensembles of trajectories at the weak excitation levels and low atomic velocities used in 
the relevant experiments. 

4. Effects of surfaces on atoms near dielectrics 

In the vicinity of a material surface, the mode structure of the full electromagnetic 
field is modified due to the dielectric properties of nearby objects. These off-resonant 
radiative interactions modify the dipole decay rate of atomic states and shift electronic 
energy levels. This surface interaction varies spatially as the relative atom-surface 
configuration changes. The surface phenomena are dispersive and depend on the multi- 
level description of the atom's electronic structure; they are calculated using traditional 
perturbation theory with the full electromagnetic field without focusing on a few select 
modes enhanced by a cavity in cQED. 

4-1- Spontaneous emission rate near a surface 

When a classical oscillating dipole is placed near a surface, its radiation pattern is 
modified by the time-lagged reflected field from the dielectric surface. The spontaneous 
emission rate oscillates with distance d from the surface, which can be interpreted as 
interference between the radiation field of the dipole and its reflection. The variation of 
the emission rate depends on whether the dipole vector is parallel or perpendicular to 
the surface, as intuitively expected from the asymmetry of image dipole orientations of 
dipoles aligned parallel and perpendicular to the surface normal. For either orientation, 
the spontaneous emission rate features a marked increase within a wavelength of the 
surface due to surface evanescent modes that become available for decay for d < Xq. 

We calculate the surface-modified dipole decay rates ^s^\d) and 'yi'^\d) for a 
cesium atom near an Si02 surface following the methods of Refs. [151 133] (see Fig. [s]). 
This calculation involves an integration of surface reflection coefficients over possible 
wavevectors of radiated light. The integrand depends on the dielectric function of Si02 
evaluated at the frequency of the atomic transition. The orientations refer to the 
alignment of a classical dipole relative to the surface plane. 

4-2. Calculation of Casimir-P older potentials 

Radiative interactions with a surface are important components of motion for neutral 
atoms within a few hundred nm of a surface, with the potential for manipulating atomic 
motion through attractive p[6] or repulsive forces [3l]. Depending on the theoretical 
framework, these forces are naturally thought of as radiative self-interactions between 
two polarizable objects, fluctuations of virtual electromagnetic excitations, or as a 
manifestation of vacuum energy of the electromagnetic field. These surface interactions, 
represented by a conservative potential f/g, are sensitive to the frequency dispersion of 
the electromagnetic response properties of the atoms and surfaces. 
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Figure 3. Variations of the dipole decay rate %{d) for a dipole oriented parallel (||) 
and perpendicular (_L) to the surface normal as a function of distance d from a semi- 
infinite region of Si02. The decay rate is in units of the vacuum decay rate 70 and the 
wavelength of the transition is A = 852 nm. 

For an atom located a short distance ci <^ Aq from a dielectric, the fluctuating dipole 
of the atom interacts with its own surface image dipole in the well-known nonretarded 
van der Waals interaction. Using only classical electrodynamics with a fluctuating 
dipole, the surface interaction potential is found to take the Lennard- Jones (LJ) form 
U^'^ = —C^/d^, where C3 is a constant that depends on the atomic polarizability and 
dielectric permittivity of the surface |35l [361 EZl EH]. At larger separations, virtual 
photons exchanged between atoms and surfaces cannot travel the distance in time 
t ~ l/u due to the flnite speed of light. Consequently, the interaction potential is 
reduced, as flrst calculated in the 1948 paper by Casimir and Polder [39] . The retarded 
surface potential takes the form f/J^* = —Ci/d'^ for a constant C4, where C4 depends 
on both c and h as this is fundamentally both a relativistic and quantum phenomenon. 
The full theory of surface forces for real materials with dispersive dielectric functions 
came with the work of Lifshitz jlOl HI]. This framework reduces to both the above 
situations for the proper limits, and, importantly, it accounts for flnite temperatures, 
predicting a Uf^ oc d^^ potential caused by thermal photons dominant at large distances 
for d ^ hc/kBT |l2]. In our discussion, we refer to these generalized dispersion forces 
as Casimir- Polder (CP) forces, whereas U^''^ , U^'^^, and U^^ refer to the appropriate 
distance limits. 

In microcavity cQED, evanescent fleld distance scales are set by the scale length of 
the evanescent fleld, Aq = Ao/27r = 136 nm (for the Cs D2 line). The relevant distances 
(0 < (i < 300 nm) span both the LJ and retarded regimes, but are much shorter than 
the thermal regime {d > 5 fim). In the transition region, the limiting power laws do not 
fully describe Us over the relevant range of d. In our modeling, we utilize a calculation 
of Us with the Lifshitz approach. The Lennard- Jones, retarded, and thermal limits arise 
naturally from the Lifshitz formalism |l2] . 
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Figure 4. Dispersive response functions for Si02 and Cesium atoms, (a) The 
dielectric function e{i^) for Si02 evaluated for frequency ^ along the imaginary axis, 
(b) Total atomic polarizability for Si02 evaluated for frequency ^ along the 

imaginary axis for the 651/2 ground state (red) and the 6P3/2 excited state calculated 
as described in [Appendix A| 



The potential Us enters our simulation in two ways. First, the transition frequency 
iOa. of the two-level atomic system shifts away from the vacuum frequency by 6^ = 
{U^^{r) — Uf{r)) / h, where Uf{f) and U^'^if) are the surface potentials for the ground and 
excited states, respectively. Since the atom transitions between the ground and excited 
states during its passage through the mode, the average net force used in calculations 
is found by weighting each contribution by the steady-state atomic state populations, 
Fs = Fi{l- (at),, + ^^(at),, (a),,. 

We calculate Uf and Uf^ for a cesium atom near a glass Si02 surface using the 
Lifshitz approach. This calculation depends on the dispersion properties of the response 
functions of materials, in this case the polarizability of the Cs ground state a{u) and 
the complex dielectric function e{uj) of the silica surface. Modeling of these functions is 
discussed in Appendix A In particular, these response functions must be evaluated on 
the imaginary frequency axis u = i^, a.s shown in Figure |4j 

Following the method of [13], curvature of the toroid surface is implemented by 
treating the toroid as a cylinder with radius of curvature R = -Dm/2. The major 
axis curvature is neglected because for all relevant distances d ^ Dm/'^- The resulting 
formula can be interpreted as a sum over discrete Matsubara frequencies ^„ = 27mkBT/ h 
with an integration over transverse wave vectors, which we quote without derivation: 



R + d' ^ 

n=0 



1 



A{R + d) 



Here, a (i^n) is the atomic polarizability and ry _|_ (i^„, k±) are the reflection coefficients of 
the dielectric material evaluated for imaginary frequency i^n- The primed summation 
implies a factor of 1/2 for the n = term. The reflection coefficients for the two 
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orthogonal light polarizations are: 



rx«„,fcx)=^^ (18) 



where 



^i + f, fcn = ^fci + e| (19) 

and e{i^n) is the complex dielectric function evaluated for imaginary frequencies i^n- 
Depending on the author, ry {r±) is sometimes referred to as ttm (^^te)- 

is calculated by numerical evaluation of (16). U^^ is also calculated using (16), 
but with an additional contribution accounting for real photon exchange from the excited 



state with the surface, which is proportional to Re 



in the LJ limit l38l 



e(^a)-l 
_e(tJa) + l_ 

The atom-surface potential for the ground state of cesium near a Si02 surface is 
shown in Fig. [5| including calculations for both a planar and a cylindrical surface. 
Without the cylindrical correction, the potential approaches the LJ, retarded, and 
thermal limits at appropriate distance scales. For the planar dielectric, our calculation 
yields C^/h = 1178 Hz /xm^ and C^/h = 158 Hz fj,m^ for the LJ and retarded limits. 
Note that the transition region between LJ and retarded regimes occurs around d ~ 100 
nm, the relevant distance scale for the experiments we are modeling. For d > D^^, the 
perturbative method accounting for the curvature is no longer accurate |13] , but at these 
distances, the surface forces are insignificant to atomic motion due to their steep power 
law fall-off. The excited state potential U^^ has a similar form to U^. 



5. Simulating atoms detected in real-time near microtoroids 

In order for the semiclassical model to be applied to our falling atom experiments, we 
must simulate the atom detection processes. In particular, in [12], falling Cs atoms are 
detected with real-time photon counting using a field programmable gate array (FPGA), 
with subsequent probe modulation triggered by atom detection. 

A microtoroidal cavity with frequency Uc is locked near the 6Si/2,F = 4 ^ 
6P3/2, -F = 5 atomic transition of Cs at uj^ at desired detuning Aca = — 00^'' ■ Fiber- 
cavity coupling is tuned to critical coupling where the bare cavity transmission vanishes, 
T < Tmin — 0.01. For atom detection, a probe field at frequency Up = Uc and flux 
Pin ~ 15 cts//is is launched in the fiber taper and the transmitted output power Pt 
is monitored by a series of single photon detectors. Photoelectric events in a running 
time window of length Atth are counted and compared to a threshold count Cth. A 
single atom disturbs the critical coupling balance so that T/Tmin > 1, resulting in a 
burst of photons which correspond to a possible trigger event. Extensive details of the 
experimental procedure are given in [12] and the associated Supplementary Information. 

Whereas only a single atom is required to produce a trigger, spectral and temporal 
data are accumulated over many thousands of trigger events since each individual atom 
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Figure 5. Atom-surface potentials (red) and U^^ (blue) for a cesium atom at 
distance d from a Si02 surface. The solid lines are for a planar surface whereas the 
dashed lines are for a curved surface with radius of curvature R — -Dm/2 = 1.5 fim. The 
limiting regimes for with a planar surface are shown as dotted lines, each calculated 
from analytic expressions not using the Lifshitz formalism. The cylindrical surface 
correction weakens the potential, which is noticeable in the retarded and thermal 
regimes. 



is only coupled to the cavity for a few microseconds. Simulation is a valuable technique 
to disentangle atomic dynamics from the aggregate data and offer insights into the 
atomic motion which underlies the experimental measurements. 

5. 1 . Simulation procedure 

Central to our simulations is the generation of a set of representative atomic 
trajectories for the experimental conditions of atoms falling past a microtoroid fulfilling 
the criteria for real-time detection. Since experimental triggering is stochastic, the 
trajectory set is generated randomly as well. For each desired collection of experimental 
parameters "P, a set of semiclassical atomic trajectories {Tj{t)}'p is generated that 
satisfies the detection trigger criteria. This ensemble is used to extract the cavity output 
functions T(t, Aap) and R{t, Aap). For each individual trajectory, t = is defined to be 
the time when the trajectory is experimentally triggered by the FPGA. For each set V, 
N is chosen large enough for a sufficient ensemble average to be obtained for the final 
output functions, which is typically at least 400 unique triggered trajectories. 

Within each simulation, the probe field is fixed to a given Wp. Cavity behavior 
is determined by the parameters cUc, k\, and Kex- h and k,\ are determined from 
measurements of the bare cavity with no atoms present. Low-bandwidth fiuctuations 
in and Uc from mechanical vibration and temperature locking are modeled as 
normally distributed random variations with standard deviations of 3 MHz and 1.5 
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MHz, respectively, that are fixed once for tlie duration of eacfi simulated trajectory. 
Similar to the experimental procedure, we impose that the bare-cavity output flux is 
less than 0.4 cts//is at critical coupling and on resonance. This rate would be identically 
zero for Acp = and critical coupling in the absence of these fluctuations. If the noise 
threshold is not met, then the particular trajectory is thrown out as it would have been 
in experiments. 

The atomic cloud is characterized by its temperature, size, and its height above 
the microtoroid. Its shape is assumed to be Gaussian in each direction with parameters 
determined by florescence imaging. For each simulation loop, an initial atomic position 
f^in is selected from the cloud and the initial velocity is selected from a Maxwell- 
Boltzmann distribution of temperature T. The trajectory is propagated forward in time 
under the influence of gravity until it crosses the toroid equatorial plane at 2; = 0. 
Only trajectories which pass within 1 fim of the toroid surface at 2; = are kept as 
a candidate for detection, as other atoms couple too weakly for triggering. Once an 
acceptable set of initial conditions is obtained, the trajectory f{t) is calculated over a 50 
fis time window around its crossing of z = 0, this time with the gravity, optical dipole 
forces, and surface interactions included. As the atom moves through the cavity mode, 
the atom-cavity coupling g, level shifts, decay rates, and forces change with position 
f, causing deviations of the trajectory from the preliminary free-fall trajectory. If the 
atom crashes into the surface of the toroid, then the coupling is set to g = onwards 
and the trajectory effectively ends (except for random 'noise' photon counts arising from 
the non-zero background transmission). 

Using f(t) and the steady-state expressions for the fields (section [2]), we find the 
transmission T(t). The photon count record Ci{tj) on each photodetector i for time 
step tj is generated from a time-dependent Poisson process with mean count per bin of 
Ci{tj) = T{tj)P[^At, where At = tj^i — tj = 1 ns and Pin is the input flux. Since the 
typical flux is Pin ~ 10 MHz and the timescale of quantum correlations is ~ 10 ns, the 
photon count process is assumed to be Poissonian on the relevant timescale of a few 
hundred nanoseconds for atom detection. The count record Ci(tj) is compared to the 
desired threshold of Cth in a time window Atth [I2]- If the trigger condition is met, the 
initial conditions rinj, Vi^j, the random cavity parameters Uc and ^ex, and the random 
number seed used to generate for diffusion processes are stored for later use. The 
semiclassical trajectory rj{t) can be fully reconstructed from these parameters. The 
time coordinate is shifted so that the trigger event occurs at t = 0. This process is 
repeated to acquire N triggered trajectories. 

Cavity output functions such as the experimentally measurable transmission 
Texp{t, V) for each simulation parameter set V are calculated from the set of trajectories 



Reflection coefficients Rexp{t,V) are calculated similarly. Spectra are calculated by 



{f,(t)}: 




(20) 
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averaging output powers over a time window ti < t < t2 for each probe frequency cUp. 
The times ti and ^2 are chosen to be the same as in our experiments, which is typically 
ti = 250 ns and t2 = 750 ns. The set of triggered trajectories {rjit)} is valid for a 
given set of conditions V and detection criteria {Cth, Atth} until the trigger at t = 0. In 
experiments, the probe frequency Up can be changed in power and detuning upon FPGA 
trigger. Although the same set of trajectories is valid before t = for each detuning, 
the trajectory set must be recalculated for t > for each probe detuning to mimic 
experimental conditions for spectral measurements. A numerical solution of the master 
equation in a number state basis is used for calculation of T{t) in (20); the linearized 
model is only used to calculate the trajectory f(t) and efficiently generate triggers. 

Whereas experiments give access only to ensemble averaged output functions, 
simulations contain the full trajectory paths. Provided that the simulation offers a 
reasonable approximation of the true ensemble of trajectories, then these results provide 
a window into the atomic dynamics underlying the cQED measurement of falling atoms 
which are not readily clear from observations. 



5.2. Simulation distributions 



The experimentally measurable cavity transmission Texp(t) is obtained in (20) as an 



average over the trajectory set {rj{t)} at each time t. Eq. (20) can formally be written as 
an integration over the probability distribution of coupling constants at time t, ptig, 0), 
for the given experimental parameters V: 

Te.p{t,V) = j dgd9T{g,9,V)Mg,9) (21) 

The function T{g, 9, V) is shown in Fig. [6]for the parameters V relevant to experiments, 
specifically with Aca/27r = 0, 60 MHz. For this discussion, we assume all frequencies 
are fixed and neglect surface shifts. In this perspective, Texp(t) is not directly related to 
the trajectory set {rj{t)} but rather the probability distribution pt{g,9) at time t. The 
time dependence of pt evolves based on the underlying trajectory ensemble. 




Figure 6. Plots of T[gtw,O^V) for (a) and (b) calculated numerically from ([t]). Atoms 
with higher ^tw generally have higher T and a larger probability for detection. The 
variation of T with 9 is evident, with a different periodicity for the two cavity detunings. 
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It is instructive to consider the probability distribution pt{g, 0) in more detail since 
it is the formal output of the simulations. We consider only the distribution Pt=o{g) 
over the coupling parameter g at the trigger time t = by integrating out the angular 



dependence. Through a reasonably simple analytic model (detailed in Appendix B), 
we calculate Pt=o{g) and compare to the results of the semiclassical simulation, which 
includes dipole and surface forces (Fig. [?]). For a cavity on resonance with the atom 
transition, Aca/27r = 0, the analytic model agrees well with a simulation when dipole and 
surface forces are not included. In this case, atom trajectories are nearly straight and 



vertical near the toroid, and the approximations of Appendix B are sufficient. When 
the full forces are included in the semiclassical model, the additional forces shift the 
distribution toward lower g. This effect is more significant for Aca/27r = 60 MHz. The 
corresponding experimental cQED spectra confirm that the semiclassical model with 
dipole and surface forces is necessary to reproduce spectral features in the real-time 
experiments (Fig. 
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Figure 7. Distributions pt=a{g) of coupling constants calculated for (a) Aca/27r = 
and (b) Aca/27r — +60 MHz. Distributions from the analytic model (red), semiclassical 
trajectory simulation with no dipole or surface forces (blue), and the simulation with 
all forces (black) are shown for comparison, (c) Experimental cQED spectra data for 
cavity detuning Aca/27r = 60 MHz (blue points) from 12J plotted with model spectra 
calculated from the distributions pt=o{g) in panel (b). The red is the analytic model 
of [Appendix B| and black is the semiclassical simulation. 



The cavity transmission T varies as a function of the atomic azimuthal coordinate 
6 = m(f), as evident in Fig. [6j This biases atomic detection towards specific locations 
around the toroid and leads to a non-uniform angular distribution pt=o{9) for atom 
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location at detection. Fig. [8] shows distributions of the atomic angular coordinate at 
the detection trigger t = for three simulation conditions relevant to the experiments 
of [12]. Although averaged spectra do not explicitly measure the coordinate 9, these 
simulation makes clear that trajectories passing through certain regions around the 
toroid are preferentially detected. The phase of the cavity output field depends on 9, 
suggesting the possibility for future experiments to measure the distribution of Fig. [8] 




9 



Figure 8. Probability distribution pt=o{0) of atomic azimuthal angle 9 ~ m(j>mod2TT 
at transit detection time t = presented as histograms of simulation runs. Shown are 
the cases for cavity detunings (a) Aca — and (b) Aca/27r = ±40 MHz. Normalization 
is such that the sum across all 6 is unity. 

5.3. Simulated trajectories 

We now turn to the simulated trajectories {rj{t)}. In contrast to experiments, in 
simulations we have the capability of turning certain forces selectively on and off. In 
particular, we can adjust the surface potential Us and the dipole forces, referred to 
symbolically as (despite them not being strictly derivable from a potential). To 
investigate the effects these optical phenomena have on atomic trajectories, we run 
simulations for four cases: the full semiclassical model, the model without surface forces 
(Us = 0), the model without dipole forces (f/d = 0), and the model without any radiative 
forces (f/d = Us = 0). 

Considering conditions relevant to [12j, we plot simulations for two sets of 
experimental parameters Vi^2 in Fig. |9} For Vi, the cavity is detuned to the red, 
whereas the cavity is blue-detuned in V2 (Aca/27r = —40 MHz for Vi and +40 MHz 
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for 7^2)- In each set of conditions, the probe field is on resonance with the cavity 
for high signal-to-noise atom detection (Acp = 0) and the average bare-cavity mode 
population of a is ~ 0.05 photons. The toroid cavity parameters are those of [T2] . 
{dmayLyh, Kin, He^}/27i = {100,11,13,17} MHz. Comparing the full model, we see that 
trajectories for Vi primarily crash into the surface, whereas those from V2 both crash 
and are repelled from the toroid. This asymmetry is due to the repulsive or attractive 
dipole force for different probe detunings relative to the atomic transition. The largest 
effect of turning surface forces off is seen in the blue-detuned trajectories, which have a 
lower crash rate when Us = 0. With f/j = 0, both Vi and V2 trajectories look nominally 
the same; the detuning A^a only affects cQED spectra, with a minor imperceptible effect 
arising from CP potentials initially shifting the atomic transition either closer to (red) 
or further from (blue) the cavity field. 

In addition to the qualitative differences in detected atom trajectories summarized 
here, the effects of and Us are evident in the experimental quantities Texp(t) and 
Rexp{t). Since here we focus specifically on trajectory calculations, the reader is referred 
to [12] for detailed comparisons of spectral and temporal simulations to experimental 
data. 
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6. Trapping atoms in the evanescent field of a microtoroid 

Our trajectory simulation can be extended to study trapping of atoms in a two-color 
evanescent far off-resonant trap (eFORT) near a microtoroidal resonator. An evanescent 
field trap takes advantage of the wavelength dependence of scale lengths for the optical 
dipole force of two optical fields with frequencies far-detuned from the atomic transition 
to limit scattering |l5l HG], HT] . The relative powers of the two fields are set so that near 
the surface, the blue-detuned, repulsive field is stronger than a red-detuned attractive 
field. As each field falls off with a decay constant of roughly A = 27r/A, at some distance, 
the red, attractive field will dominate and the atom will be attracted to the surface 
forming a potential minimum. Recently, evanescent fields have been harnessed to trap 
atoms in a two-color eFORT around a tapered optical fiber [IS] , where the fiber enables 
efiicient optical access to deliver both high intensity trapping fields and weaker probe 
fields to the trapped atoms in a single structure. The tapered fiber can be positioned 
as desired, bringing the trapped atoms near a device for atomic coupling. 

The tapered nanofiber eFORT is a remarkable achievement toward integrating 
atom traps with solid-state resonators, but the nanofiber scheme does not allow direct 
integration with a cavity for achieving strong, coherent coupling between light and 
trapped atoms. Another disadvantage is that trap depth is limited by the large total 
power required to achieve trapping with evanescent fields. The high quality factors 
and monolithic structure of WGM resonators allow evanescent field traps free from 
these problems while maintaining efiicient optical access from tapered fiber coupling. 
Two-color evanescent field traps in WGM resonators have been analyzed in detail for 
spheres [IH] and microdisks [50J. In this section, we extend our simulations of atoms 
in the evanescent field of a microtoroid to an eFORT that can capture single falling Cs 
atoms triggered upon an atom detection event. 

Unlike nanofibers, a microtoroid cannot be placed directly in a magneto-optical 
trap for a source of cold atoms. As shown in [12], we have the experimental capability 
to detect a single atom falling by a microtoroid and trigger optical fields while that atom 
remains coupled to the cavity mode. The semiclassical simulations described here are 
ideal for investigating the capture of falling atoms in a trap triggered upon experimental 
atom detection. 

We add an additional eFORT potential Ut to our semiclassical trajectory model in 
addition to the dipole forces and surface potential f/g. For our simulation, Ut is formed 
from a red (blue)-detuned mode near 898 nm (848 nm) with powers ~ 50 fiW to give a 
trap depth of ~ 1.5 mK at c? ~ 150 nm from the surface (Fig. [lO^). The red (blue) fields 
interact primarily with the 66*1/2 — ^ QP1/2 {GS1/2 — )■ 6P3/2) transition. The trap depth 
is limited by the total power in vacuum that can propagate in the tapered fibers of [12] . 
Power handling can be improved with specific attention to taper cleanliness, so with 
experimental care the trap depth can be increased reasonably from the discussion here, 
although we simulate under the conditions given to illustrate that this trap is already 
experimentally accessible. 
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Figure 10. (a) The trapping potential Ut along the z = axis with the CP potential 
included. Also shown are the red and blue evanescent potentials of the two trapping 
modes, Ut, respectively, (b) The mode function used in Ut for the 898 nm mode 
with m = 106. (c) Simulated trajectories for trapping simulations with a eFORT Ut 
triggered "on" by atom detection at i = with Aca = 0. Falling atoms with the FORT 
beams "off' {t < 0) are colored blue, whereas trajectories after the trap is triggered 
are red. Trajectories are colored pink for t > 50 /is to illustrate the timescale. Roughly 
25% of the triggered trajectories become trapped, (d) Same as (c) showing only the 
trapped trajectories and a clearer view of atom orbits in the evanescent trap. 



The difference in vertical scale lengths {ipo in Q) for modes of different wavelength 
leads to a trap that is not fully confined if both the red- and blue-detuned trap modes 
are of the lowest order (as in Fig. [2l1o). As increases, the repulsive blue-detuned 
light weakens faster than the red-detuned field, and atoms can crash into the toroid 
surface. This problem is alleviated by exciting a higher-order mode for the 898 nm 
light, as shown in Fig. \T0\ ). The modal pattern confines atoms near z = and prevents 
trap leakage along tp. This problem is not present in the microdisk eFORT of [50] 
because the optical mode extent is determined by structural confinement and not the 
optical scale length. Use of a higher-order mode was also used to form an atom-gallery 
in a microsphere j49j . 

During the detection phase of the simulation, Ut = 0. At t = conditioned on an 
atom detection trigger, Ut is turned on. The kinetic energy of an atom with typical fall 



Simulations of atomic trajectories near a dielectric surface 



24 



velocity of i; ~ 0.2 m/s is equivalent to 0.3 mK, so a 1.5 mK trap is sufficiently deep to 
capture an atom if it is triggered near the trap potential minimum. Defining a trapped 
trajectory to be one such that the atom has g/2n > 5 MHz at t = 10 /is, approximately 
25% of triggered atom trajectories are captured when the trapping potential is turned 
on. Simulated trapping times exceed 50 fis, limited not by heating from trapping light 
but by the radiation pressure from the unbalanced traveling whispering-gallery modes 
of nearly-resonant optical probe field. This probe field can be turned off so that the 
atoms remain trapped beyond the simulation time. 

In contrast to the standing-wave structure of a typical eFORT or Fabry-Perot cavity 
trap [5T], microtoroidal resonators offer the tantalizing possibility of radially confining 
an atom in a circular orbit around the toroid [HI [52]. The Ut = outlined here 
does not confine the atoms azimuthally, forming circular atom-gallery orbits around the 
microtoroid [52] (Fig. 10 :,d) . In the same manner as jlH] , a localized trap can be achieved 
by exciting a red-detuned standing wave for three-dimensional trap confinement. 

This trapping simulation outlines how real-time atom detection can be utilized to 
trap a falling atom in a microtoroidal eFORT. In practice, microtoroidal traps present 
some serious practical challenges. Notably, because the trap quality is sensitive to the 
particular whispering-gallery mode, the excited optical mode must be experimentally 
controlled. The success of an eFORT for Cs atoms around a tapered nanofiber |1H] 
strongly suggests that similar trap performance might be achieved for an eFORT 
around a high-(5 WGM cavity, localizing atoms in a region of strong coupling to a 
microresonator. 



7. Conclusion 



We have presented simulations of atomic motion near a dielectric surface in the regime of 
strong coupling to a cavity with weak atomic excitation. As required by experimental 
distance scales, this simulation includes surface interactions, which manifest through 
transition level shifts and center-of-mass Casimir-Polder forces. Analysis of the 
simulated trajectories gives insight into the atomic motion underlying experimental 
measurements of ensemble-averaged spectral and temporal measurements for single 
atoms detected in real-time. We have adapted our simulations to investigate the 
capturing of atoms in an evanescent field far off-resonant optical trap in a microtoroid. 
Our simulations suggest that falling atoms can be captured into an eFORT around 
a microtoroid, offering an experimental route towards trapping a single atom in atom- 
chip trap in a regime with both strong cQED interactions and significant Casimir-Polder 
forces simultaneously. In this system, the sensitivity afforded by coherent cQED can be 
used not only for atom-chip devices, but also as a tool for precision measurements of 
optical phenomena near surfaces. 
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Appendix A. Calculating the Polarizability and Dielectric Response 
Functions 



Evaluation of Casimir-Polder interactions of atoms with the surface of the dielectric 
resonator requires evaluation of the atomic polarizability and of the dielectric function 
as functions of a complex frequency. Here we outline our analytic model of the complex 
dielectric function for Si02 and the atomic polarizability of Cesium atoms in the ground 
and excited states. 

The complex dielectric function e^u) = ei+ie2 is modeled using a Lorentz oscillator 
model of the real and imaginary parts of the response function to analytically introduce 
frequency dependence and enforce causality, 

e{u) = e^ + J2 1^2^ 2 2 (K' - + ^^7,) (A.l) 

J J J 

Here, Uj is the resonance frequency, 7^ is the damping coefficient, and fj is the oscillator 
strength for each oscillator in the model, e^o = e(a; — )■ 00) = 1. e can be expressed in 
terms of the complex index of refraction n = + m as e = = — + 2inK, where 
n is the refractive index and k is the extinction coefficient. Experimental data for h for 
Si02 is available over a wide frequency range which is used to fit the parameters 



of (A.l) for a seven-oscillator model {j = 1 — 7). Using the analytic form of (A.l), the 
dielectric function can readily be evaluated over complex frequencies as shown in Fig. |4j 
The frequency-dependent atomic polarizability as{ijj) for Cesium in a state s is 
calculated as a sum over transitions of the form, 

where e is the electron charge, me is the electron mass, Uns is the transition frequency, 
and fns is the signed oscillator strength for the transition of state n to the state s 
ifns > if state n is above s in energy). A more complete expression for the response 
function a{u) should include damping coefficients given by the transition linewidths. 
Since our calculations involve integrals over infinite frequency on the imaginary axis 
and atomic linewidths are generally narrow with respect to transition frequencies, we 
assume that the off-resonant form given by (A. 2) without damping is sufficient. We also 
note that this expression does not account for the differences between magnetic sublevels 
and hyperfine splitting, which again represent small corrections when these expressions 
are integrated over the imaginary frequency axis. The general form of (A. 2) applies to 
the polarizabilities for both the 6S1/2 ground state and the 6P3/2 excited state, with an 
additional tensor polarizability for the QP3/2 state. 
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The total atomic polarizability is composed of contributions from valence electron 
transitions (a^) and high-energy electron transitions from the core shells to the 
continuum (etc), such that a = + oic- The valence polarizability constitutes 96% 
of the total static polarizability (5^ in Cs, with only significant at high frequencies. 
We take Oc to be the same for both the ground and excited states of Cs, whereas a„ 
is obviously sensitive to the different electronic transition manifolds for 65*1/2 and 6P3/2 
states. Valence electron oscillator strengths and transition frequencies are tabulated 
in many sources |55l [56]. Our estimate of a^{uj) for the ground state includes all 
65'i/2 — > iVPi/2 and 6S'i/2 — > ^^-^3/2 transitions, with = 6 — 11. For the excited 
state, a^iuS) is calculated using 6P3/2 — )■ (6 — 15)S'i/2, 6P3/2 — t- (5 — 11)1^3/2, and 
6P3/2 ~^ (5 — ll)-D5/2 transitions. Tensor polarizability contributions sum to zero 
when averaged over all angular momentum sublevels [HZ]. In agreement with [31], our 
calculation of ctv comprises about 95% of the total static polarizability. 

For simplicity, all core electron transitions are lumped into a single high-frequency 
term of the form used in (A. 2). This term contains two free parameters, /core and cucore, 
which are found from the following two conditions. Using the calculation of a^iiS) for 
the Cs ground state, we enforce that the ground state static polarizability aiuj — )■ 0) 
matches the known value calculated theoretically [58] a(0) = 5.942 x 10~^^ cm'^. We 
also ensure that the ground state LJ constant for a Cs atom near a metallic surface 
agrees with the known value [Ml [59] C3 = - ^ OL{ii)di = 4.4 ■ h kHz /im . These 
conditions are sufficient to fix the two free parameters in ac{u}) for this single oscillator 
core model, although the high-frequency structure of the core polarizability is lost. For 
the excited state calculation, we use the same ar(co)- 



Appendix B. Analytic model of falling atom detection distributions 

Here we develop an analytic model of the distribution p{a.n{g, of coupling parameters 
g and azimuthal coordinate 6 = mcf). Atoms are assumed to fall at constant vertical 
velocity with no forces, in contrast to the more complete semiclassical trajectories used 
in this manuscript to generate Pt{g)- An abbreviated description of this model appears 
in the Supplementary Information of |12j . 

The linearized steady-state cavity transmission T(Aap, is a known function of 
Aap and r. We only consider the lowest order mode where the cavity mode function 
is approximately Gaussian in z and exponential in distance from the surface d. The 
approximate temporal behavior of the coupling constant g for a single trajectory is, 

^7(p,^(t))=^c(p)e-(^«/^°)^ (B.l) 

where gdp) is the maximum value of the g at the closest approach of its trajectory 
(2; = 0), zq is a characteristic width assumed to be independent of p, and z{t) = —vt. 
gdp) decays exponentially from the maximum (yfmax at the toroid surface, gdp) ~ 
fi'maxe^^'^^^P''^^". The transmission T and hence the detection probability depend on 
6] in general, if atoms fall uniformly around the toroid, the most numerous trajectories 
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detected will be at the values of 6 which maximize T{6) for the cavity parameters of 
interest {9 = tc/2 for Aca/27r = +40 MHz, for example, as in Fig. |8]). 

The probability density function for the full ensemble of detected falling atoms 
Piaiiig, 9) can be estimated as the product of the probability of any atom having a 
particular g and the probability of a trigger event occurring for an atom with coupling 

9, 

Pi^\\{.9, 9) ~ Patom(fi')Ptrigger(5', &) . (B.2) 

An atom transit is triggered when the total detected photon counts exceeds a threshold 
number, Cth, within a detection time window Atth- For a probe beam of input flux 
Pin, the mean counts in this window are C = T{g, 6)PinAtth- This expression assumes 
that the atom is moving slowly so that the T{g, 9) at trigger event is the only T{g, 9) 
that contributes to the detection probability. The detection probabihty Ptrigger(5', 6*) is 
estimated from a Poisson distribution of mean count C . 



From (B.l), Patom(fi') can be written as a product of the probability p{g\gc) of an 
atom in a trajectory with a given g^ to have coupling g and the probability of a trajectory 
to have that g^, Pmaxigc), integrated over all gc, 

Patom(fi') = / P{.g\gc)Prn^^{gc) dgc (B.3) 

•'9 

The integral has limits from g to (^max since gc cannot be smaller than g. 

For atoms falling uniformly over the p — 4> plane, Pmax(5'c) dgc is proportional 
to the area of a ring of radius p and thickness dp, Pma^igc) dgc ~ 2Tipdp. Using 
gc{p) ~ e-^P-^'-l^^l^^ ^ ~ - J. Hence, pmax(^?c) ~ l/gc for (p - Dp/2) < Z^p/2. 
To find p{g\gc) we note that that the probability is proportional to the time an atom in 



the trajectory is at a particular g. From (B.l) for a constant velocity f, this trajectory 



is Gaussian and the relative probability must be proportional to dz. Finding the 
differential as a function of g gives p{g\gc) oc dz 



Putting the results together in (B.3) gives 



9\/^K9c/9) 



P.toM - I — ^f= - V^^^— (B.4) 

ggc ^Hgjg) g 



9 

This result diverges as g goes to zero since there are infinite transits with small gc and 
infinite time for atoms with small g for any transit regardless of gc for t — )■ ±oo. This 
divergence is not problematic in calculating (B.2) since Ptrigger(5'' ^) '^^^^ ^^'^ 9 
faster than the logarithmic divergence in Patoml^ ■ 

The spectrum for given experimental parameters as a function of probe detuning 
Aap = Up — wi^'' can be written as: 

r(Aap) = / r(Aap, g, 9)pi^n{g, 9) dg d9 (B.5) 

where the normalization of Pfsiii{g,9) is chosen such that 



5'niax 



Pmi9,9)dgd9 = l (B.6) 
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The overall probability of g, Pfa.\i{g) independent of 9, is found by integrating over 9. 
In practice, P{a.n{g) is quite similar to Pfa.ii{g,9) evaluated for the 9 which maximizes 
the transmission. Fig. [t] compares this simple model for pisiiig) with the equivalent 
distribution from the semiclassical trajectory simulation, pt=o{g). 
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